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a^ ■ 

0^ ' We introduce a model of Poincare mappings which represents hierarchical structure of 

^-H ^ phase spaces for systems with many degrees of freedom. The model yields residence time 

< I , distributions of a power type, and hence temporal correlation remains long. The power law 

' behavior is enhanced as the system size increases. 

OO ' Introduction In many Hamiltonian systems, l/f^^O < v < 2) power spectra ^ 

■ long time tails have been observed, for instance, in^rea preserving mappings, 
I a water cluster, and a ferro-magnetic spin system. The 1/ spectra imply that 

J> ' relaxation to equilibrium is slow. They are hence important phenomena of Hamilto- 

. nian systems with many degrees of freedom. We are interested in understanding the 

Q I cause of 1//^ spectra from the structure of phase space and properties of motion. 

' To describe 1/ spectra, Aizawa introduced a geometrical model for area pre- 

serving mappings, which are models of Poincare mappings for Hamiltonian systems 
with two degrees of freedom. u^'EP This model assumes exact self-similar hierarchical 
structure of phase spaces and produces stagnant motion, namely slow relaxation. 
J>^. We therefore understand that stagnant motion arises from self-similar structure of 

T3 I phase space (often referred to ed as "islands around islands" Q)) and motion tranped 

Q ■ to KAM tori or Cantori. Meiss et al. have successfully proposed a similar model. 13 A 

^ . renormalization group approach, B*'© which demonstrates similarity between scale 

I transformations in ph&sg space and in time also supports the picture described above. 

• However, the models H^' EP are based on the two-dimensionality of the phase space and 

. cannot be directly applied to high dimensional systems. 

^ I For systems with many (more than two) degrees of freedom, Aizawa et al. 

■ discussed the origin of the 1/ spectra based on the Nekhoroshev theorem. Since 
the argument is based on the Nekhoroshev theorem, the relation between stagnant 
motion and the hierarchical structure of phase space is not clear. 

Moreover, the assumptions on which the models mentioned above are based do 
not seem to hold for high-dimensional systems. There exist many sorts of fixed 
points of Poincare mappings from the fully elliptic type to the fully hyperbolic type, 
and only the fully elliptic fixed points yield exact self-similarity as area preserving 
mappings. Since the ratio of fixed points of the fully elliptic type decreases as the 
system size becomes large, it seems impossible to assume exact self-similarity in 
the phase space structure for general high-dimensional Hamiltonian systems. It is 
believed that KAM tori rapidly disappear as the systems size becomes large, and 
hence the volume of the region where stagnant motion occurs also decreases. On the 
other hand, since stagnant motion is frequently observed for Hamiltonian systems 
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with many degrees of freedom, we need to establish a model which yields stagnant 
motion in systems with many degrees of freedom accordingly. 

In this paper, we propose a geometrical model which represents hierarchical 
structure of phase spaces. The model is an extension of Aizawa's model to many de- 
grees of freedom, and assumes that sticky zones exist around fixed points of Poincare 
mappings even if fixed points are not fully elliptic. In other words, in systems with 
N degrees of freedom, motion is assumed to be trapped for a time around tori of 
fewer than N dimensions also. 

Types of fixed points We consider Poincare mappings F and their fixed points 
instead of Hamiltonian flows and their periodic orbits. We set the number of degrees 
of freedom to n for Poincare mappings which have 2N — 2 dimensional Poincare 
sections, where N is degrees of freedom of the Hamiltonian dynamics (i.e. n = N—1). 
Note that we construct a model based on fixed points of F hereafter. We can 
construct the model based on periodic points of F with period-/c by using F'^ instead 
of F. 

Local structure around fixed points is built from a combination of the three 
elementary types: elliptic, hyperbolic and vortex types, for which the eigenvalues 
of Jacobian of F are (e**^, e"*"!), (r, 1/r) and (re*"^, re~*'^, e*'^/r, e~*'^/r) respectively, 
where both u> and r are real. lilP The local structure is constructed as direct products 
of these three types. 

We assume that there is no vortex type structure for simplicity. Generalization 
including the vortex type will be given in Ref. [T^ ). Then local structure around a 
fixed point is constructed by elliptic and hyperbolic types of structure, and there are 
re + 1 varieties of fixed points: direct products of re — i elliptic type and i hyperbolic 
(z = 0, 1, • • • ,n). We define the index of a fixed point as i. For instance, a fully 
elliptic fixed point is index-0 and a fully hyperbolic is index-re. 

Geometrical model and master equation Let us introduce a geometrical model with 
the following assumptions: 

(G-1) Hierarchical structure is constructed by fixed points of Poincare mappings in 
phase spaces. 

(G-2) Every sort of fixed points has a sticky zone around it, even if it is not fully 
elliptic. 

We calculate volumes of each level of hierarchy and derive a master equation with 
some assumptions. The number of the level is put in order by volume, and the base 
level is level-0. We assume that the regions of level-(/ -|- 1), {I + 2), - ■ ■ are in the 
region of level-/. The schematic picture of this model is described in Fig. 1. 
Let us introduce notation for the quantities which will be used later: 

Ni^i : number of fixed points of level-Z, index-i, 

Vi^i : volume of a sticky zone around a fixed point of level-/, index-i, 

Ni : total number of fixed points of level-/, 

Vi : total volume of sticky zones of level-/, 
pi^i : number of sets of fixed points of level-(/ -|- 1), 
which surround a fixed point of level-/ index-i. 
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Fig. 1. Schematic picture of the hierarchical structure of phase space. In this picture we assume the 
system size to be n = 2. Black points are fixed points. Circle and cross under the fixed points 
represent elliptic and hyperbolic elementary types, respectively. The double circle, for instance, 
implies that the structure around the fixed point is a direct product of two elliptic types, namely 
index-0. Let us focus on the left half of this figure. The fixed point at the center, whose index 
is 0, belongs to level-Z, and the other four fixed points belong to level-(Z + 1). Shaded areas 
are sticky zones around fixed points of level-(Z + 1). The sticky zone of the fixed point at the 
center is the inside of the biggest circle, which includes sticky zones of level-(i + 1) fixed points. 
Level-(Z + 2) fixed points exist around level-(Z + 1) fixed points, with a similar situation as on 
the right half, which is magnification around a fixed point with level-(Z + 1) and index-0. 

Here we have assumed: 

(G-3) Fixed points have the same volume Vi^i if their level and index are the same. 
Note Vi > Vi+i for all / because we require that Vi includes V;+i, V;+2, • ' The 
meaning of pi^i is clarified in the following. 

(Fact): A fixed point of level-/ and index-i is surrounded by pi^i{n-iCj) fixed 
points of level-(/ + 1) and index- (z where j = 0, 1, • • • , n — i, and pi^i is a positive 
integer. 

This fact is an extension of the Poincare-Birkhoff theorem El-* for many degrees 
of freedom. EJ) Using this fact, we write the recursion formula for A^; j as 

i 

Nl+l^i = pl^kNl^k{n-kCi-k)- (1) 
fc=0 

The total number Ni and volume Vi of the level-Z sticky zones are 

n n 

Ni = J2n^, Vi = Y.Ni^,Vi,i. (2) 

i=0 1=0 

To observe motion among levels we introduce a master equation with the fol- 
lowing three assumptions: 
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(M-1) Systems are ergodic. 

(M-2) Transitions from level-/ are limited to level-(/ — 1), / and (/ + 1). 
(M-3) Transitions among levels are Markovian. 

From (M-1), the probability being level-/ in equilibrium, Pi"^, is proportional to 
volume V/: 

pr yi- (3) 

We assume detailed balance to fix the transition probability Wm,i = w{l m}. 
Then 

wi^^P^ = w^jPp. (4) 
From (M-2) and Eqs. (^ and (Q), transition probabilities are written as 

cVi^i, (m = / - 1) 

cVi+i, [111 = 1 + 1) 

0. (otherwise) 

The factor c is independent of the level, and we set c = l/(2(Vo -|- V2)). By using 
these transition probabilities, the master equation is written as 

Pl{t + l)=Y,Wl,mPmit), (6) 

m 

where Pi{t) is the probability of being on level-/ at step t. 

Results of numerical calculations We numerically calculated residence time dis- 
tributions based on our model Eq. (^) and examined if it obeys a power law. The 
residence time distribution R(t) is the probability that motion extends to levels shal- 
lower than level-/th for the first time with initial level being /th at i = 0. We obtain 
R{t) from the following transition probabilities and initial condition: 

Wm,i^^-i = (m = /th - l,/th,/th + 1), Pi{t = 0) = 61^1^^. (7) 

Here ^1,1^^^^ is Kronecker's delta. Then we have 

Rit) = Pi^^-i{t). (8) 

If R{t) is of a power type rather than an exponential, motion among levels is stagnant. 
Parameters which must be given are assumed as follows: 

a, = 2, Vi,i = b-^'+'\ /th = l, 

iVo,o = 1 and iVo,i = 0. (i > 0) 

Sticky zones become small as hyperbolic components increase, and this index de- 
pendence of Vi^i is an essential point of our model. We assumed that the volume of 
a sticky zone Vi^i decays as b~'^ with respect to the index-i of the fixed point, since 
local structure near fixed points is constructed by direct products. Other forms of 
Vi^i, for instance Vi^i = b^''b~^ and Vi^i = exp[— (/ + i°')\nb], give similar results to 
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Fig. 2. Residence time distributions for various values of the scale factor b. Here n = 80. The 
magnitude of the longitudinal axis is multiplied by 10"\ 10"^ and 10"^ for 6 = 40, 50 and 60, 
respectively. 

V;^j = 6~('+*) with appropriate values of parameters. The form of VJ^j determines 
which index is dominant in V;. 

The residence time distribution -R(t) is shown in Fig. 2 for n = 80 and various 
values of h. When h is small (h = 30, 40) is of a power type. Namely, we have 

R{t) ~ t-^ (9) 

where /3 is 3.4 and 3.1 for 6 = 30 and h = 40, respectively. Since R{t) is of a 
power type, stagnant motion occurs among levels. The stagnant motion is also 
observed when n = 1, and hence our model Eq. (^) is consistent with models for 
area preserving mappings mentioned in the Introduction. The values of /3 are also 
consistent with those obtained in area, pre ser ving mappings which are models of 
physical systems and give 1.5 < /3 < 3.I1P'd''I13) As far as we know, /? has not been 
calculated in Hamiltonian systems with many degrees of freedom. Appearance of 
oscillating behavior for 6 = 50 and 6 = 60 is caused by the weakness of effects 
of the hierarchy, which become weaker as b increases, because large b implies that 
level-(^ + 1) is small compared with level-L 

We display the residence time distribution R{t) for various system size n in Fig. 
3, where power law behavior of R{t) is clearly seen. Oscillations found for small n 
gradually decay as n becomes large, and the distributions are close to a power type 
as n increases. This is an effect of the many degrees of freedom and indicates that 
fine tuning of parameters is not necessary to observe l/f^ spectra in systems with 
many degrees of freedom. 

Summary To understand l/f^ spectra and long time tails in Hamiltonian systems 
with many degrees of freedom, we proposed a geometrical model of phase space, 
which is an extension of Aizawa's model. We assumed that sticky zones exist around 
fixed points of Poincare mappings even if the fixed points are not fully elliptic, and 
accordingly, exact self-similarity of phase space is not introduced. We derived a 
master equation from our model, and found that residence time distributions are of 
a power type. That is, stagnant motion among levels occurs although phase space 
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Fig. 3. Dependence on the number of degrees of freedom of residence time distributions with fixed 
b = 30. The magnitude of the longitudinal axis is multiplied by lO'S 10"^, 10"^ and 10"* for 
n = 20, 40, 60 and 80, respectively. The distributions are closer to a power type as n increases. 

does not possess exact self-similarity. The power law behavior becomes clearer as 

the system size increases. 

We express our thanks to members of R-lab. of Nagoya University for useful 
discussions. 
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